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ABSTRACT 


GRB 221009A (z = 0.151) is one of the closest known long y-ray bursts (GRBs). Its extreme 
brightness across all electromagnetic wavelengths provides an unprecedented opportunity to study a 
member of this still-mysterious class of transients in exquisite detail. We present multi-wavelength 
observations of this extraordinary event, spanning 15 orders of magnitude in photon energy from radio 
to 4-rays. We find that the data can be partially explained by a forward shock (FS) from a highly- 
collimated relativistic jet interacting with a low-density wind-like medium. The jet’s beaming-corrected 
kinetic energy (Ex ~ 4 x 10°° erg) is typical for the GRB population, but its opening angle (~ 2°) is 
one of the narrowest. The radio and mm data provide strong limiting constraints on the FS model, but 
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require the presence of an additional emission component. From equipartition arguments, we find that 
the radio emission is likely produced by a small amount of mass (€ 6 x 1077 M) moving relativistically 
(T > 9) with a large kinetic energy (> 104° erg). However, the temporal evolution of this component 
does not follow prescriptions for synchrotron radiation from a single power-law distribution of electrons 
(e.g. in a reverse shock or two-component jet), or a thermal electron population, perhaps suggesting 
that one of the standard assumptions of afterglow theory is violated. GRB 2210094 will likely remain 
detectable with radio telescopes for years to come, providing a valuable opportunity to track the full 


lifecycle of a powerful relativistic jet. 


1. INTRODUCTION 


Long-duration 7-ray bursts (GRBs) release enormous 
amounts of energy in the form of collimated, highly- 
relativistic jets. These outflows are thought to be 
launched by a powerful magnetar or accreting black hole 
born during the collapse of a massive star (Piran 2005; 
Woosley & Bloom 2006; Marshall et al. 2011). GRBs are 
typically discovered via their prompt (seconds — minutes 
duration) y-ray emission, possibly powered by internal 
shocks within the jet (Rees & Meszaros 1992; Kobayashi 
et al. 1997; Kumar & Piran 2000). Properties of the 
jet, such as its structure, energetics, and magnetization, 
along with other physically interesting quantities such 
as the structure and density of the ambient medium, 
can be derived from modeling the broadband afterglow 
emission, which typically remains detectable in the ra- 
dio, optical, and X-ray bands for days to months (e.g., 
Sari et al. 1998). 

Synchrotron models, based on both analytical and 
numerical or hydrodynamic calculations, have been re- 
markably successful in explaining the multi-wavelength 
afterglow emission of GRBs. In the basic picture, the 
afterglow is modeled as synchrotron radiation produced 
by à population of relativistic electrons accelerated in 
the forward shock (FS) produced by the interaction of 
the jet with the ambient medium. The radiating elec- 
trons are assumed to be accelerated into a single power- 
law distribution of energies, characterized by a minimum 
energy, YmMec? and a power-law index p. This creates a 
simple broken power-law emission spectrum, fully char- 
acterized by a few break frequencies and overall flux 
normalization. Physical model parameters such as the 
jet energy and opening angle, the microphysical shock 
parameters, and the density and density profile of the 
circumburst medium may be determined by capturing 
the full synchrotron SED and its temporal evolution, 
which typically requires observations spanning the full 
electromagnetic spectrum. 

Despite the success of this basic picture, increasingly 
detailed observational campaigns of GRBs over the past 
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decade have demonstrated the need to incorporate ad- 
ditional physics. One of the most important additions 
has been the detection of reverse shock (RS) emission, a 
second synchrotron component from a second shockwave 
propagating back into the jet (Laskar et al. 2013; Perley 
et al. 2014; van der Horst et al. 2014). The RS emission 
reveals details of the jet's magnetization and initial bulk 
Lorentz factor (Sari & Piran 1999; Kobayashi & Zhang 
2003; Zhang & Kobayashi 2005). Additional model ex- 
tensions, from the inclusion of multiple ejecta shells 
moving at different Lorentz factors (Rees & Meszaros 
1998; Jóhannesson et al. 2006; Laskar et al. 2015) to 
the consideration of emission from electrons not accel- 
erated in the FS (thermal electrons; Eichler & Waxman 
2005; Ressler & Laskar 2017; Warren et al. 2018; Mar- 
galit & Quataert 2021), have also been explored in the 
literature. Further extensions may be needed to eluci- 
date several aspects of GRB observations that remain 
unexplained. These include the detection of very high- 
energy photons (> 1 TeV) associated with some GRBs 
(Abdalla et al. 2019; MAGIC Collab. 2019; Blanch et al. 
2020; H. E. S. S. Collaboration et al. 2021) and unusual 
radio evolution inconsistent with variability due to ex- 
trinsic scattering from the interstellar medium or the 
standard synchrotron emission framework (Frail et al. 
2004; Bright et al. 2019; Kangas & Fruchter 2021; Levine 
et al. 2023). Numerical simulations may provide the ul- 
timate solution to these oddities, but rely on uncertain 
assumptions about particle acceleration physics, mag- 
netic field structure, etc. 

On 2022 October 9, the Fermi Gamma-Ray Burst 
Monitor (GBM) and the Swift Burst Alert Telescope 
(BAT) triggered on a new y-ray source, rapidly con- 
firmed to be the brightest GRB ever seen by these in- 
struments and designated GRB 221009A (Dichiara et al. 
2022; Veres et al. 2022; Williams et al. 2023). These de- 
tections were followed by numerous others across the 
electromagnetic spectrum, confirming a bright optical, 
X-ray, and radio counterpart (Dichiara et al. 2022; 
Lipunov et al. 2022; Kennea et al. 2022a; Bright et al. 
2022) and a redshift of 0.151 (de Ugarte Postigo et al. 
2022; Malesani et al., in prep.). GRB 221009A also ex- 
hibited coincident very high-energy emission, with pho- 


tons reported up to 18 TeV (Huang et al. 2022) and 
potentially as high as 250 TeV (Dzhappuev et al. 2022). 
Although a calculation of the isotropic-equivalent energy 
released in y-rays for this burst is complicated by satu- 
ration effects at multiple y-ray instruments, preliminary 
estimates yield log Ey iso ~ 54.77 (0.1 keV-100 MeV), 
one of the highest values to date (Kann & Agui Fernan- 
dez 2022). 

Several groups have modeled the afterglow emission 
from GRB 221009A (Ren et al. 2022; Sato et al. 2022; 
Guarini et al. 2023; Fulton et al. 2023). Fulton et al. 
(2023) focus solely on the optical and X-ray afterglow 
in the context of characterizing a possible emerging su- 
pernova (SN) component at ~ 20 days. The other three 
groups use radio, optical, X-ray, and/or y-ray observa- 
tions taken at < 12 days post-discovery (Ren et al. 2022; 
Sato et al. 2022; Guarini et al. 2023). Although the de- 
tails of each model differ, common findings among the 
three groups include a high isotropic-equivalent kinetic 
energy and a high degree of collimation for the jet and 
a low-density environment. In all analyses, the best-fit 
model cannot explain all of the data, suggesting that 
GRB 221009A may be an ideal test case for more com- 
plex models incorporating more realistic physics. 

Here, we present detailed radio, millimeter, optical, X- 
ray, and y-ray observations of GRB 221009A, spanning 
fifteen orders of magnitude in frequency and four orders 
of magnitude in time (extending to ~ 100 d). We find 
that a standard synchrotron FS model can explain large 
portions of the data, revealing that GRB 221009A’s ex- 
traordinary brightness is due to the jet’s unusual degree 
of collimation, rather than an intrinsically large energy. 
However, the model struggles to reproduce some aspects 
of the full dataset, in particular the radio and millimeter 
emission. We show that standard extensions to the syn- 
chrotron FS model (e.g., the addition of a reverse shock) 
are unsuccessful at explaining our observations, suggest- 
ing that at least one of the basic assumptions underlying 
our standard picture of GRB afterglow emission needs 
to be revisited. 

Unless otherwise noted, all times refer to times af- 
ter the GBM trigger (2022 Oct 09 13:16:59.000 UTC; 
Veres et al. 2022), and all magnitudes are in the AB 
system, not corrected for Galactic or intrinsic dust ex- 
tinction. We employ ACDM cosmological parameters 
of Qm = 0.31, Q4 = 0.69, and Ho = 68 kms~! Mpc^! 
throughout, and assume a Galactic extinction of Ay = 
4.1034 mag (Schlafly & Finkbeiner 2011). The redshift 
of z = 0.151 corresponds to a luminosity distance of 
Dry, = 2.28 x 102” cm for this burst. 


2. OBSERVATIONS AND DATA ANALYSIS 


2.1. Radio 
2.1.1. GMRT 


We observed the afterglow using the upgraded Giant 
Metrewave Radio Telescope (GMRT) through program 
43.039 (PI: Laskar) starting 18.0 days, 21.9 days, and 
21.8 days after the burst in Bands 3 (400 MHz), 4 
(750 MHz), and 5 (1260 MHz), respectively. The 
observations utilized 400 MHz bandwidth. We used 
J1925+2106 for complex gain calibration, and either 
3C 286 or 3C 48 for bandpass and flux density cali- 
bration, depending on availability. For Band 3, we car- 
ried out data reduction via standard techniques using 
the Common Astronomy Software Applications (CASA; 
McMullin et al. 2007), including three rounds of phase 
self-calibration, followed by two rounds of amplitude and 
phase self-calibration. For Band 4 and Band 5, we re- 
duced and imaged the data using the capture pipeline 
(Kale & Ishwara-Chandra 2021). We list the results of 
our uGMRT observations in Table 1. 


2.1.2. MeerKAT 


We observed the afterglow with the MeerKAT radio 
telescope beginning at ~ 1.3 days after the burst in 
L Band (1.3 GHz) via program SCI-20220822-TL-01 
(PI: Laskar). Later observations employed simultaneous 
subarray-mode observations in the L and UHF bands. 
The observations employed 1938-634 as bandpass and 
flux density calibrator, and J1925+2106 as complex gain 
calibrator. The pipeline SDP images were of sufficient 
quality for photometry of the afterglow using CASA. 
The noise in the image of L-band data taken on 2022 
Oct 15 was much higher than in the other epochs due to 
significant RFI contamination. We report our MeerK AT 
photometry in Table 1. 


2.1.8. VLA 


We obtained seven epochs of multi-frequency observa- 
tions of GRB 221009A beginning ~ 3.5 days after the 
burst with the Karl G. Jansky Very Large Array (VLA) 
through program 22B_062 (PI: Laskar). The observa- 
tions employed either 3C 286 or 3C 48 for bandpass 
and flux density calibration and J19254-2106 for com- 
plex gain calibration. We imaged the pipeline-calibrated 
measurement sets downloaded from the VLA archive, 
where available. For the second epoch taken on 2022 
Oct 15, the flux calibrator scans were defective. We 
calibrated this observation using flux calibrator data at 
the same frequencies taken from the first epoch instead, 
and achieved excellent agreement for the derived flux 
density of the gain calibrator. The pattern of the imag- 
ing residuals in the second and third epochs indicated 
gain errors, which we ameliorated via phase and ampli- 
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tude self-calibration using an automated self-calibration 
pipeline!. We report the results of our VLA observa- 
tions in Table 1. 


2.2. VLBA 


We obtained four epochs of observations of 
GRB 221009A with the Very Long Baseline Array 
(VLBA) and the 100-m radio telescope Effelsberg at 
8.3 GHz beginning ~ 12.4 days after the burst un- 
der project VLBA/22B-305 (Legacy Code BLO73; PI: 
Laskar). The full-track observations utilized 3C 345 
and 3C 454.3 as fringe calibrators and were phase refer- 
enced to TXS 1903+196 (JVAS J1905+1943), a compact 
quasar at z = 2.3, with cycles of 168s on target and 
68s on the calibrator. We carried out standard fringe 
fitting?, bandpass, and complex gain calibration using 
AIPS (Greisen 2003). We interpolated the delay and 
rate solutions to GRB 221009A and produced phase- 
reference images. We measured the flux density using 
JMFIT in AIPS in the image plane and with MODELFIT in 
DIFMAP (Shepherd et al. 1994) in the uv plane, obtain- 
ing similar results within a few percent. We corrected 
the derived flux densities for the primary beam for the 
first epoch and re-scaled according to the calibrator 
flux in the fourth epoch. We present the inferred total 
(CLEANed) flux density measurements from our VLBA 
observations in Table 1. 


2.2.1. ATCA 


We obtained four epochs of observations of the afterglow 
with the Australia Telescope Compact Array (ATCA) 
beginning z 5.8 days after the burst by triggering our 
program C3289 (PI: Laskar) using the CABB correla- 
tor tuned to 15 mm (with 2 GHz basebands tuned to 
16.7 and 21.2 GHz) and 7 mm (33 and 35 GHz; and 
also 45 and 47 GHz). The observations were carried 
out by the observatory in service mode under project 
CX515, and employed PKS1921-293 for bandpass cali- 
bration, PKS1934-638 for flux density calibration, and 
PKS1923+210 (J1925--2106) for complex gain calibra- 
tion. We calibrated the data using MIRIAD (Sault et al. 
1995), and combined the 33 GHz and 35 GHz data, as 
well as the 45 GHz and 47 GHz data, prior to imaging 
in CASA. The cleaning process reveals residuals charac- 
teristic of phase decorrelation and the resulting ATCA 
SEDs are extremely steep and are inconsistent with an 


1 https://github.com/jjtobin/auto_selfcal 


2 To account for the structure of the complex gain calibrator, we 
derived a clean-component model by first imaging the data af- 
ter a preliminary fringe fit and used this model for subsequent 
calibrations. 


extrapolation to nearly contemporaneous ALMA obser- 
vations. Whereas phase-only self-calibration recovers 
some flux, the signal-to-noise in the data is too low 
for adequate self-calibration. We list the ATCA mea- 
surements in Table 1 for completeness, but we do not 
use these in our subsequent modeling, and we caution 
against the use of these data in other works without 
more careful attention to the calibration. 


2.2.2. ALMA 


We obtained seven epochs of ALMA Band 3 (3 mm) 
observations of GRB 221009A beginning ~ 2.4 days 
after the burst through program 2022.1.01433.T (PI: 
Laskar). The observations utilized two 4 GHz wide base- 
bands centered at 91.5 and 103.5 GHz, respectively. The 
first five epochs used J1924-2914 as flux density calibra- 
tor. The sixth epoch used J1550+0527; this execution 
was affected by a correlator issue and was re-observed, 
this time using J22324-1143. The seventh epoch used 
J1550+0527. All epochs used J1914+1636 as complex 
gain calibrator. Our ALMA coverage has a gap be- 
tween z 15-75 days due to the shutdown of the ob- 
servatory following a cyber attack. We downloaded the 
pipeline-generated images from the ALMA archive and 
performed photometry in CASA. We report our pho- 
tometry in Table 1. 


2.2.3. NOEMA 


We obtained two epochs of NOEMA 97.5 GHz obser- 
vations of GRB 221009A at 39.2 and 54.2 days after 
the burst through program S22BE (PI: Laskar). The 
observations utilized two 7.7 GHz wide basebands cen- 
tered at 89.8 GHz and 105.2 GHz, respectively. We used 
MWC2349 as flux density calibrator, B2200+420 (epoch 
1) and 3C 454.3 (epoch 2) as bandpass calibrators, and 
1932+204 and 1923+210 as complex gain calibrators. 
Data were reduced in GILDAS? using standard proce- 
dures by observatory staff and provided to us. We im- 
aged the reduced data in CASA and report our photom- 
etry in Table 1. 


2.2.4. SMA 


We observed the afterglow with the Submillimeter Ar- 
ray (SMA) at a combination of 1.3 mm (~ 230 GHz) and 
1.1 mm (290 GHz) for 7 epochs using Uranus as flux 
density calibrator, 3C 84 as bandpass calibrator, and 
interleaved observations of J1925+211 and MWC349a 
for complex gain calibration. We calibrated the data in 
MIR‘ (the in-house calibration suite for the SMA) and 


3 https://www.iram.fr/IRAMFR/GILDAS 
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'Table 1. Radio observations of GRB 221009A 


Telescope Frequency? Time? Flux density Uncertainty 
(GHz) (days) (mJy) (uJy) 
MeerK AT 1.28 1.26 2.10 24 
ALMA 97.5 2.36 9.21 90 
SMA 231 2.68 9.44 900 


NOTE—? Central frequency. P Mid-time since Fermi/GBM trigger. 


The full data table is available on request. 


measured the afterglow flux density using vector aver- 
aging of the uv data (verified by imaging of the first two 
epochs). An additional 850 um observation yielded an 
upper limit. We report the results of our SMA observa- 
tions in Table 1. 


2.3. Ultraviolet 


We performed photometry in a 5” aperture on all 
Swift/UVOT images up to and including segment 
01126853067 (taken on 2022 Dec 16) obtained from the 
Swift repository? with the uvotproduct (v2.8) software 
and CALDB version 20221229. This corresponds to al- 
most all of the data taken before the target entered a 
Sun constraint on 2022 Dec 21 (lasting until 2023 Feb 
06). We used defaults for all pipeline parameters. We 
present our UVOT measurements in Table 2 (see also 
Williams et al. 2023 for an independent analysis of these 
data). The white-band data is of limited utility owing 
to the strong foreground extinction and we therefore do 
not use data in this band for subsequent modeling. 


Table 3. Liverpool Telescope Observations of 
GRB 221009A 


AL? Filter Mag Uncertainty Tit Seeing 


(d) (AB) (s) (arcsec) 
0.33486 g 18.57 0.04 45 1.17 
0.33749 Z 15.29 0.01 45 1.04 
0.33662 i 15.96 0.01 45 1.08 
0.33576 i 16.99 0.01 45 1.40 
0.34899 g 18.53 0.04 50 1.01 
0.34987 r 17.05 0.01 40 0.99 


'Table 3 continued 
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Table 2. 
GRB 221009A 


Swift/UVOT Observations of 


Start Time? Stop Time? Band Mag Uncertainty 


(s) (s) 
593.2 612.9 UVB 17.09 0.13 
767.0 786.7 UVB 17.21 0.14 
1147.9 1332.8 UVB 17.50 0.13 
40736.7 41644.4 UVB 20.21 0.14 
56651.5 57558.5 UVB 20.72 0.26 
337.2 586.9 UVU 17.67 0.08 
742.1 761.9 UVU 17.78 0.25 
1122.9 1317.9 UVU 17.77 0.18 
28831.1 29661.9 UVU 20.3 0.21 
46079.7 46907.8 UVU 20.61 0.27 
92151.0 92979.8 UVU 20.83 0.31 
668.4 688.3 UVV 15.56 0.11 
841.2 861.0 UVV 15.55 0.11 
1049.4 1243.1 UVV 15.80 0.09 
34595.2 35419.7 UVV 18.07 0.07 
51699.8 52523.8 UVV 18.64 0.10 
80607.5 81514.3 UVV 19.10 0.13 
97761.4 98526.1 UVV 19.54 0.20 


NoTE—?Since Fermi/GBM trigger. Magnitudes are in the 
native Swift/UVOT system and have not been corrected for 
extinction in the Milky Way or GRB host galaxy. 


Table 3 (continued) 


AL? Filter Mag Uncertainty tint Seeing 
(d) (AB) (s) (arcsec) 
0.35061 i 16.00 0.0 30 0.90 
0.35132 Zz 15.33 0.0 30 0.89 
0.39745 g 18.55 0.05 50 .36 
0.39978 z 15.48 0.0 30 .02 
0.39834 r 17.21 0.0 40 .01 
0.39908 i 16.15 0.0 30 0.96 
1.27230 r 18.85 0.02 45 0.87 
1.27097 g 20.33 0.06 100 .34 
1.27404 Z 17.06 0.0 45 .03 
1.27316 i 17.77 0.0 45 0.84 
1.30908 r 18.87 0.02 90 .08 
1.31046 i 17.80 0.0 90 .03 
1.31186 Zz 17.10 0.0 90 -00 
1.39474 g 20.41 0.12 100 iL? 
1.39607 r 19.03 0.04 45 .06 
1.39694 i 17.93 0.0 45 .41 
1.39781 Zz 17.21 0.0 45 14 
2.26839 i 18.69 0.02 100 ll 
2.26689 r 19.84 0.05 100 .35 
2.26991 Zz 18.02 0.0 100 -02 
2.29593 r 19.84 0.05 60 04 
2.29441 g 21.05 0.12 120 .33 
2.29698 i 18.73 0.02 60 1.24 


Table 3 continued 


Table 3 (continued) 


At * Filter Mag Uncertainty tint Seeing 


(d) (AB) (s) (arcsec) 
2.29803 z 7.99 0.02 60 1.24 
3.27665 i 9.29 0.03 120 1.51 
3.27434 r 20.31 0.05 180 1.58 
3.27863 Z 8.43 0.07 120 2.69 
3.28182 g 21.61 0.19 180 2.18 
3.28371 r 20.33 0.09 60 1.91 
3.28475 i 9.28 0.03 60 1.46 
3.28579 Zz 8.55 0.03 60 1.56 
7.27887 id 21.60 0.08 180 0.79 
7.28154 i 20.44 0.03 180 0.85 
7.28571 Z 9.68 0.03 360 0.75 
7.28992 i 20.37 0.03 180 0.83 
7.29263 r 21.44 0.08 180 0.82 
8.27806 i 20.58 0.04 240 1.07 
8.27470 r 21.68 0.08 240 1.07 
8.28328 Zz 19.94 0.04 480 1.20 
8.28854 i 20.69 0.05 240 1.08 
8.29195 r 21.91 0.13 240 1.09 
10.28995 r 22.04 0.11 300 0.90 
10.29366 i 21.07 0.06 240 0.84 
10.29887 z 20.27 0.05 480 1.04 
10.30412 i 20.92 0.05 240 0.85 
10.30788 È 22.04 0.10 300 0.86 
11.27403 r 22.09 0.11 360 0.83 
11.28344 Z 20.50 0.05 480 0.80 
11.27821 i 21.12 0.05 240 0.83 
11.28871 i 21.17 0.06 240 0.91 
11.29292 r 22.09 0.10 360 0.91 
14.28998 r 22.45 0.17 360 0.90 
14.29496 i 21.47 0.07 360 0.87 
14.30262 Z 20.83 0.06 720 0.79 
14.31111 i 21.35 0.09 240 0.95 
14.31535 r 22.35 0.15 360 0.99 
17.31003 Z 21.08 0.15 720 1.67 
18.29124 i 21.91 0.13 720 1.37 
18.30665 r 22.99 0.24 1350 1.46 
19.27566 Zz 21.19 0.73 600 1.08 
20.27930 z 21.55 0.14 1350 1.23 
21.28883 z 21.31 0.12 1350 1.23 
36.28571 i 22.23 0.11 1620 1.10 
38.27 709 Ë 24.06 0.32 1800 0.87 


NoTE—?Mid-time since Fermi/GBM trigger. The data have not been 
corrected for extinction in the Milky Way or GRB host galaxy, or 
for the contribution of host galaxy light. 


2.4. Optical 


We observed GRB 2210094 in griz filters with IO:O 
on the Liverpool Telescope (LT; Steele et al. 2004) at 
multiple epochs beginning on 2022 Oct 09. We down- 
loaded pipeline-reduced images from the LT archive. To 
avoid contamination from nearby sources in the crowded 
field, we use a custom script to construct a model of the 


PSF and subtract the wings of several bright stars after 
masking the inner pixels, as needed. We report aperture 
photometry performed in a 1.0" radius aperture cali- 
brated to the Pan-STARRSI catalog in Table 3. We ad- 
ditionally include optical and near-infrared (NIR) data 
reported in GCN circulars, from Xu et al. (2022); Brivio 
et al. (2022); Rastinejad & Fong (2022); O'Connor et al. 
(2022a); Liu et al. (2022); Huber et al. (2022); Ferro 
et al. (2022); O'Connor et al. (2022b) and radio obser- 
vations from Leung et al. (2022) in our analysis. We do 
not include Hubble Space Telescope (HST) observations 
from Levan et al. (2022) in our work as the correspond- 
ing photometry likely contains contamination from an 
underlying host galaxy, and require deep, late-time tem- 
plates for subtraction. 


2.5. X-rays and y-rays 
2.5.1. NuSTAR (3-79 keV) 


The Nuclear Spectroscopic Telescope Array (NuS- 
TAR, Harrison et al. 2013) acquired a first set of four 
observations between 1.96 — 23.73 days (PIs Margutti 
& Racusin) with exposure times of ~ 20ks, followed 
by one deeper exposure at 31.91 days (exposure time of 
40 ks, PI Troja). We reduced the NuSTAR data follow- 
ing standard procedures with the NuSTAR Data Analy- 
sis Software (NuSTARDAS) version 0.4.9 and NuSTAR 
CALDB version 20221130, applying standard filtering 
criteria with nupipeline. GRB 221009A is a bright 
source of hard X-rays in the NuSTAR bandpass (3-79 
keV) at all times. For each epoch we extracted a spec- 
trum with nuproducts using a source extraction region 
centered at the location of the radio counterpart and dif- 
ferent sizes to maximize the S/N as reported in Table 4. 
We used a source-free background region of radius > 2'. 
We find that the hard X-ray spectrum is well modeled 
by a simple power-law with photon index? of r ~ 1.86 
at all times. 


2.5.2. Swift/XRT (0.3-10 keV) 


The Swift X-ray Telescope (XRT) began observing 
GRB 221009A z 3.4 ks after the Swift/BAT trigger and 
& 6.6 ks after the Fermi/GBM trigger. We extracted 
XRT PC-mode spectra at the times corresponding to the 
NuSTAR epochs using the time-sliced spectrum tool on 
the Swift website’ and modeled the spectra, together 
with corresponding calibration files, in XSPECv12.12.1. 


6 We define the photon index with the convention, Fg x E- T. 
7 https://www.swift.ac.uk/xrt_spectra/01126853/ 


Table 4. Log of NuSTAR observations 


Start Date/Time Stop Date/Time 


Centroid MJD . Time? 


Net Exposure? Net Exposure? Source region size 


(dd-mm-yy/hh:mm:ss) | (dd-mm-yy/hh:mm:ss) (d) 
2022-10-11/10:04:09 2022-10-11/14:45:00 1.96 
2022-10-15/05:21:09 2022-10-15/17:16:09 5.92 
2022-10-20/01:06:09 2022-10-20/11:56:09 10.72 
2022-11-02/06:06:09 2022-11-02/17:01:09 23.93 
2022-11-09/23:06:09 2022-11-10/23:01:09 31.91 


(d) A (ks) B (ks) Radius (7) 
1.96 20.66 20.49 1.5 
5.92 20.66 20.49 1.1 
10.72 20.44 20.26 1.0 
23.73 21.30 21.09 1.0 
31.91 40.78 40.38 1.0 


NoTE—* With respect to the Fermi/GBM trigger time. 
> For NuSTAR module A. 
* For NuSTAR module B. 


While the derived photon index appears to increase with 
time (from z 1.55 to z 1.82, tying Nyint across epochs), 
the evidence for this is marginal (S 2c) and we do 
not consider this statistically significant (but see also 
Williams et al. 2023). Tying the photon index across 
epochs gives results consistent with the parameters on 
the Swift website. 

Finally, we combine the XRT and NuSTAR data to- 
gether and perform three joint spectral fits: (i) tying 
both Ny int and Dx across epochs; (ii) tying Nuint and 
allowing [x to vary; and (iii) tying Lx and allowing 
Nu,int to vary. We do not find evidence for evolution 
in [x with time. There is marginal evidence in these 
fits for a decrease in Ny in, from 1.51792 x 102? cm? 
to 0.99*03 x 102? cm-?. However, these numbers are 
consistent at the 2c level, and we do not consider vary- 
ing Ng, further. On tying both quantities across 
epochs, we find Nyint = 135006 x 102? cm7? and 
Dy = L8566 0055. 

We download the XRT 0.3-10 keV count rate light 
curve ? (Evans et al. 2007, 2009) and convert it to flux 
density at 1 keV for subsequent analysis (after shift- 
ing the time to the GBM trigger time) using [xy = 
1.8566. For this we use unabsorbed counts-to-flux con- 
version rates of 1.00 x 1071? ergcm"? ct^! and 1.07 x 
1071? erg cm? ct^! for the WT and PC-mode, respec- 
tively, obtained from the Swift website. We also extract 
the NuSTAR flux in the range 15-20 keV and convert it 
to flux density at 15 keV using the photon index from 
the joint fit for subsequent analysis. 


2.5.8. Fermi/LAT (100 MeV-100 GeV) 


The LAT instrument on board the Fermi satellite is 
sensitive to gamma-ray photons in the energy band from 
30 MeV to 300 GeV (Atwood et al. 2009). We extracted 


8 Obtained from the Swift website at http://www.swift.ac.uk/xrt_ 
curves/01126853 and rebinned to a minimum signal-to-noise ratio 
per bin of 10. 


and analyzed the Fermi/LAT data of GRB 2210094 us- 
ing the public software GTBURST, which is distributed 
as part of the official Fermitools software package?. We 
extracted the LAT data within a temporal window from 
3.5 ks to 100 ks (0.04-1.16 days) after the GBM trigger 
time. We filtered photons with energies in the 100 MeV— 
100 GeV range, within a region of interest (ROI) of 12? 
centred on the burst position of R.A. = 288.264? and 
Dec. = 19.773°, and with an angle from the spacecraft 
zenith < 100°, as part of the standard procedure. We 
selected the “P8R3 SOURCE” class as the instrument 
response function, suitable for late-time emission subse- 
quent to the prompt phase of the burst. We extracted 
light curves assuming a power-law model for the spec- 
trum of the source, together with the “isotropic tem- 
plate” and “template” for the particle background and 
the Galactic component, respectively. We performed an 
unbinned likelihood analysis, setting the minimum test 
statistics (TS) to 10 for detections. We calculated flux 
upper limits using a spectral slope of —2 in the energy 
band of 100 MeV-100 GeV. We report our Fermi /LAT 
photometry in Table 5. 


3. MULTI-WAVELENGTH MODELING 


We now interpret the multi-wavelength afterglow ob- 
servations in the context of the standard model of syn- 
chrotron radiation produced by a relativistic forward 
shock (FS) produced by the GRB jet propagating into 
the pre-explosion environment. The model is parame- 
terized by the shock energy, Ex iso, the radial density 
profile p x r-^ of the environment, the density nor- 
malization (ng for k — 0, ISM-like; and A, for k — 2, 


9 https:/ /fermi.gsfc.nasa.gov /ssc/data/analysis/scitools/gtburst. 
html 


Table 5. Fermi/LAT Observations 


Start Time? Stop Time? Energy Flux Energy Flux Photon Flux Photon Flux Photon Index’ Photon Index Test 
Uncertainty Uncertainty Uncertainty Statistic 
(s) (s) (ergs^! cm?)  (ergs^! cm?)  (phs^!cm^?)  (phs-^!cm-?) 

3000.0 4429.2 2.67 x 10-05 1.43 x 10-95 3.06 x 10-95 1.06 x 107° 2.10 0.283 57 
4429.2 6539.4 9.54 x 10-9? 2.69 x 10-9? 1.12 x 10-95 1.88 x 10-06 2.11 0.141 177 
9654.9 14254.6 4.44 x 10-99 1.73 x 10-9? 5.41 x 10-96 1.37 x 10-98 2.13 0.209 79 
14254.6 21045.8 1.72 x 10-9? 9.79 x 10-19 2.35 x 10-96 9.06 x 10-97 2.19 0.334 30 
21045.8 31072.3 1.44 x 10-99 7.95 x 10-19 2.01 x 10-98 8.81 x 10-97 2.20 0.354 21 
31072.3 45875.7 5.70 x 10-19 2.79 x 10-19 id" 10-99 5.26 x 10-97 2.46 0.402 13 
45875.7 67731.6 4.34 x 10-1? 1.80 x 10-1? 1.27 x 107° 5.08 x 10-97 2.87 0.445 12 
67731.6 100000. < 6.42 x 107 1° ET < 5.79 x 10-07 s 2.00 5 


NoTE—* With respect to the Fermi/GBM trigger time. > Here defined as Fg x ETT. 
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Figure 1. Light curves of GRB 221009A from 1 GeV to 400 MHz together with the broken power-law model fits used to 
interpolate data to a common time and to inform the multi-wavelength modeling (Sections 3.1 and 3.2). The light curves have 
been multiplied by the factors listed in the legends to offset them in the plot for clarity. The three panels with radio data are 
on the same temporal scale to facilitate comparison. 


wind-like!?), the fraction of shock energy in relativis- 
tic electrons (€e) and magnetic fields (ep), as well as 
the index (p) of electrons accelerated to a power-law en- 
ergy distribution. The observed radiation is expected 
to be characterized by three break frequencies: (i) the 
self-absorption frequency, v4; (ii) the characteristic fre- 
quency, Vm; and (ii) the cooling break, ve. Solving 
for these parameters requires interrogating the observed 
light curves and spectral energy distributions (SEDs). 


3.1. X-ray and optical/NIR: Forward Shock 


In order to interpolate observations to common times 
for constructing and plotting SEDs, we fit the observed 
multi-wavelength light curves to a series of broken power 
law models. The resulting fits are presented in Fig- 
ure 1. The Fermi/LAT 1 GeV and NuSTAR 15 keV 
light curves can be fit with a single power law with decay 
indices!!, auar = —1.4740.08 and ays key = 1.744:0.03. 
The Swift/XRT 1 keV light curve can be fit with a bro- 
ken power law with decay indices, ax, = —1.46 + 0.02 
and ox,» = —1.70+0.03; however, the break time, ty, x = 
(1.00.8) day is not well constrained. The r'-band light 
curve can be fit with a broken power law with decay in- 
dices, a,,; = —0.48 + 0.06 and a, 2 = —1.40 + 0.02 and 
break time, tp, = (0.32 + 0.04) days, although the pre- 
break decay is contingent upon uncertain photometry 
reported in GCN circulars. Correcting the r'-band flux 
for Milky Way extinction, the NIR-to-optical spectral 
index at ~ 4.4 days is BNIR—opt = —0.76 + 0.04, while 
the NIR-to-X-ray spectral index is similar, ÜNiR x = 
—0.70 + 0.01. Spectral indices involving NIR/optical 
bands are subject to additional, indeterminate system- 
atic errors due to unknown uncertainties on the Galactic 
extinction along the line of sight and any additional in- 
trinsic extinction. 

A shallower light curve in the optical compared to 
the X-rays usually indicates!? an ISM-like environment 
with Vopt < vc < vx. However, a uniform-density envi- 
ronment is ruled out for this burst by the closure rela- 
tions between the light curve decay rates and spectral 
indices. To see this, we start with the optical spectral 
index of ÊNIR-opt ~ —0.76, which requires an electron 
energy index of p z 2.52 in the regime Vopt < vc < vx. 
This would imply ox — (2 — 3p)/4 z —1.4, which is 
shallower than observed. Furthermore, the similarity 
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between NIR-opt and nrr-x suggests that no spectral 
break is present between these bands. If we instead con- 
sider the regime opt < vx < Ve, and p z 2.52 from the 
spectral indices, then the predicted light curve decline 
rates of Qopt © ax & —1.1 are significantly shallower 
than observed in either band. Finally, in the regime 
Ve < Vopt < vx, the observed spectral indices would 
indicate p ~ 1.5, which predicts a light curve decline 
rate of a = —(3p + 10)/16 ~ —0.9 (Dai & Cheng 2001; 
Zhang & Mészáros 2004), again much shallower than ob- 
served in either band. While this can be remedied by 
interposing a jet break (after which the light curve de- 
clines at & —p, including exponential sideways spread- 
ing; Rhoads 1999; Sari et al. 1999) at tjet 7: 1 day, in 
this model the jet quickly becomes non-relativistic and 
the resulting evolution cannot match the radio obser- 
vations (Section 3.2). We note that Sato et al. (2022) 
came to a similar conclusion even without radio data, 
but incorporating Fermi/LAT GeV observations. For 
completeness, we present an ISM model, together with 
its shortcomings, in Appendix A. 

We next consider the other standard case of k — 2 
for the density profile of the pre-explosion environment. 
The ve < Vopt < vx regime is ruled out as in the ISM 
case, as in this regime the light curves are agnostic to 
the density profile. In the regime Vopt < Ve < vx, the 
optical spectral index again requires p ~ 2.52. The pre- 
dicted light curve decay rates are Qop 7 —1.64 and 
ax & —1.39; however, the opposite is observed!?. In 
addition, in this regime we expect an X-ray spectral 
index of Bx = —p/2 ~ —1.26, which is significantly 
steeper than the observed. Finally, we consider the 
regime Vopt < vx « vc. The spectral index between the 
Swift/XRT band at 1 keV and the Fermi/LAT band 
at 1 GeV is Byey_-cey = —1.10 + 0.04. This is sig- 
nificantly steeper than the X-ray spectral index alone. 
Furthermore, the spectral index in the GeV band of 
Baev = 1 — Faev © —1.2 (using a weighted-mean pho- 
ton index, Table 5) also appears steeper than the X-ray 
data alone, supporting the presence of a spectral break 
between the X-ray and GeV bands. While the optical 
and X-ray spectral indices and decay rates ought to be 
equal to each other in this regime, a potential solution! 
for the shallower optical decay observed could be the 
proximity of Vm to the optical/NIR bands. 


10 4, = 1 corresponds to a mass-loss rate of My = 107? Mo/yr for 
a wind velocity of vw = 10? km s 1. 

11 We use the convention F, œ t*v? throughout. 

1? Whereas the presence of a SN could also make optical light curves 
shallower, as we discuss later in this section, we find no strong 
additional evidence for such a component. 


13 This remains a fundamental issue in modeling this burst, as we 
discuss in Section 4.3. 

14 Fulton et al. (2023) assume the same spectral regime of Vopt < 
Vx < ve and instead explain the shallower optical light curve 
compared to the X-rays as contribution from an underlying SN 
in the optical. However, they also require strongly time-varying 
extinction for this scenario to explain the observed optical color. 
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SEDs (left) and a sampling of available light curves (right) at X-ray, UV, optical, and NIR wavelengths of 


GRB 221009A (data points), together with the highest-likelihood FS wind model (lines) that explains most of these obser- 
vations (Section 3.1). Thin lines indicate the underlying synchrotron spectrum. The break at ~ 2 x 10'* Hz in the first epoch 
is/$. The X-ray SEDs shown here are the result of joint spectral fits of Swift and NuSTAR observations at the corresponding 
times (Section 2.5.2). Optical points have been interpolated using the best-fit broken power-law light curves (Figure 1). The 
model cannot explain the UV/optical light curves at < 1 day; these points (open symbols) were excluded from the fit. These 


caveats are discussed further in Section 4.3. 


We find that we are able to find a relatively satis- 
factory solution with k = 2 for p ~ 2.5 in the regime 
Vm Í Vopt < Vx < Vc < VpyaT (Figure 2) with the follow- 
ing caveats: (i) this model does not match any of the 
radio observations (although this is also true for the ISM 
model, see Appendix A), and over-predicts the ALMA 
97.5 GHz data point at ~ 100 days; (ii) the model over- 
predicts the NIR (JHK) flux at the ~ 30% level; and 
(ii) in the optical, this model under-predicts the ob- 
served light curves at both < 1 day and = 30 days. We 
discuss point (i) in Section 3.2 and points (ii) and (iii) 
in Section 4.3. 

To explore this FS model further, we fit the 
Fermi/LAT light curve at 1 GeV, the Swift/XRT light 
curve at 1 keV, the UV /optical/NIR data at 1-30 days, 
and the 97.5 GHz data point at ~ 100 days simultane- 
ously by sampling the parameter space of p, E iso, Ax, 
€e, €B, and tje using Markov Chain Monte Carlo with 
emcee (Foreman-Mackey et al. 2013). The details of 
our implementation are described in Laskar et al. (2013, 
2014). We run 512 walkers for 5k steps and discard the 
first 100 steps as burn-in. We use a uniform prior on p 
from 2.001 to 2.99. We restrict log(e.) and log(eg) to the 
range € (—5,0) and require e, + ep < 1. We constrain 
log(A,) € (710,10), Ex;iso/erg € (1048,5 x 1055), and 
log(tjet/days) € (—5,5) and use Jeffreys (1946) priors 
for these last five parameters. To account for systematic 
flux calibration offsets in any given band in our data set 
as well as to prevent single, high signal-to-noise points 
from driving the entire fit, we implement a minimum un- 


certainty floor of 1096 prior to running the fit. We fix the 
Galactic extinction to Ay ga] = 4.1 mag, and find some 
evidence for additional extinction, Ay host % 0.2 mag. 
However, we note that there is a degeneracy between 
Ay,gal and Ay host as the redshift is low; thus it is en- 
tirely possible that the true Galactic extinction is lower 
and the true extinction along the line of sight through 
the host galaxy is higher than the inferred value. 

In the highest-likelihood model, the synchrotron break 
frequencies are located at va zz 2 x 10° Hz, Vm ~ 4.5 x 
101^ Hz, and ve ~ 3.8 x 10!? Hz at 1 day, with Vm c 
Vopt < Vx < Vc < VpAT, as expected. The proximity of 
Vm to the optical contributes to a shallower optical light 
curve than expected for the regime Vm < Vopt < wc. 
The proximity of v, to the X-rays results in a spectral 
index intermediate between £8 = (1 — p)/2 ~ —0.76 and 
B = —p/2 ~ —1.3. The resulting SED fits the NuSTAR 
spectra well, even though those data were not included 
in the fit. Finally, the spectral index of 8 = —p/2 = 
— 1.3 closely corresponds to the value of fGey % —1.2 in 
the Fermi/LAT band at 1 GeV. 

We summarize the results of our MCMC analysis in 
Table 6. The electron index p is sharply constrained by 
the X-ray light curve, the NIR-to-X-ray spectral index, 
and the X-ray-to-GeV spectral index. As v4 is uncon- 
strained in this model, the physical parameters ( Ex iso; 
Ax, €e, €B) exhibit degeneracies (Figure 3). The ob- 
served 3mm flux at ~ 100 days is lower than expected 


from a spherical FS, and this drives a jet break!? at 
around this time in the model, yielding a small jet open- 
ing angle, Ojet S 2°. 

As mentioned earlier, this model (i) somewhat (~ 
30%) over-predicts the NIR (JH K) flux (Figure 2) and 
(ii) under-predicts the optical emission at S; 1 day and 
at Z 30 days. We speculate that the mismatch in the 
NIR may be related to the inability of this model to also 
fit the radio data, as discussed next in Section 3.2. We 
discuss these discrepancies in the NIR and optical, to- 
gether with additional caveats on the FS modeling, in 
Section 4.3. 

We note that our model parameters are somewhat dif- 
ferent from the analyses of Ren et al. (2022) and Sato 
et al. (2022) (notably, we find a much higher value of 
cpg). However, these papers rely solely on data collected 
within S 7 and < 12 d of GRB 221009A’s discovery, 
respectively. When compared to our more extensive 
dataset, we find that both previously proposed models 
dramatically overpredict the radio emission at t 2 5 d 
and are therefore ruled out. Finally, in contrast to Ful- 
ton et al. (2023), we do not include a SN contribution 
to our model, as this does not appear to be required by 
the data. In particular, there is no strong evidence for 
excess emission relative to the afterglow model in any 
of the optical or NIR light curves. This suggests that 
the contribution of the SN is lower than the flux level of 
the observed multi-frequency optical/NIR light curves. 
Further investigation of the SN requires spectroscopic 
information, and we defer a detailed discussion of the 
SN to papers focused on this emission component. We 
also do not include host galaxy emission in our model. 
While host emission is known to affect the HST data 
(not included here) at ~ 30 days (Levan et al. 2022), we 
do not expect the host galaxy to make a significant con- 
tribution in the optical light curves at the earlier times 
(S 20 days) considered here. 


3.2. Radio: Multiple Components 


In Figure 4, we present our VLA SEDs together with 
data from GMRT, ALMA, and NOEMA at 400 MHz, 
800 MHz, 1.3 GHz, 97.5 GHz, and 225 GHz interpo- 
lated! to the times of the VLA SEDs using their corre- 
sponding best-fit broken power-law functions (Figure 1). 
We also extract XRT SEDs at the times of the VLA 
SEDs and fit for the spectral normalization with spec- 


15 We note that the jet break time is bounded below by the absence 


of a steep decline in the X-ray light curve to tjet Z 70 days. 


16 The typical temporal dispersion of data points in a given SED is 


~ 1096. 
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Table 6. FS Model Parameters 


Parameter Best fit MCMCÓ 

p 2.531 2.53 + 0.01 
log €e —0.616 | —0.56*007 
log ep —0.221 — —0.40*0 26 
log A, -3.012  —3.50*0 19 
log( Ex iso /erg) 54.04 53.99 + 0.07 
log(tjet/d) 2.110 2.141019 
Viet 1.46 1.647935 
Av, host 0.24 0.23+ 0.02 
log( Ex /erg) 50.55 50.60* 002 
log va t /H2)? 5.3 

log(va,r/Hz) 14.6 

log vc,c/ Hz) 19.6 

log( F3, /mJy) 0.98 


NoTE—Frequencies and flux densities are cal- 
culated at 1 day. 


"Summary statistics from the marginalized 
posterior density distributions, with median 
and +34.1% quantiles (corresponding to 

+1o for Gaussian distributions; Figure 3). 


b This frequency is not directly constrained by 
the data. 


tral parameters fixed from the NuSTAR-XRT joint fit 
(Section 2.5.2). 

We find that the FS model discussed in Section 3.1 
under-predicts all the radio observations, except at the 
highest frequencies (at 97.5 GHz, from ALMA) at = 
75 days. The radio emission is strongly self-absorbed 
below + 2 GHz for most of the period spanning ~% 3- 
76 days, while the optically thin spectrum above the ra- 
dio peak does not match the inverted v!/? spectrum ex- 
pected from the low-energy tail of the minimum energy 
electrons. In addition, the spectral index above the peak 
(B ~ —0.2) is significantly shallower than GNtR—opt as 
well as 8x. This indicates that either the radio emission 
arises from a separate emission component, or the ap- 
proximation of the electron power-law being truncated 
at a minimum Lorentz factor (ymin) breaks down in this 
case. 

Fitting the radio SEDs at S 50 GHz with fiducial, 
broken power-law models of the form 


/s 


—1 
1 v —sßpı " 1 V —sf»2 
2 Vpeak 2 Vpeak 


(1) 
separately in each epoch, we find evidence for a slowly 
decreasing peak flux (fy peak o t-0 7950.02) and peak 
frequency (Vpeak cx t7 9495902) with time (Figure 5 and 


F, = Lf peak 
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Figure 3. Correlations and marginalized posterior density from multi-wavelength modeling of GRB 221009A with Ex,so and 
Ex in erg. We compute the opening angle (jet) in degrees from Ex iso, Ax, and tjet- The contours enclose 39.3%, 86.4% and 
98.9% of the probability mass in each correlation plot (corresponding to lo, 2c, and 3e regions for two-dimensional Gaussian 
distributions), while the dashed lines in the histograms indicate the 15.996, 5096 (median), and 84.196 quantiles (corresponding 
to +1lo for one-dimensional Gaussian distributions). The distribution for eg cuts off at high values due to the prior bound of 
€e + ep < 1. Strong correlations are evident between most parameters as v4 is unconstrained in this model. See Table 6 for a 


summary. 


Table 7). From these fits, we confirm that the spec- 
tral index above the peak is shallow, Bradio ^: —0.2 
at < 28 days, steepening marginally to ~% —0.4 at 
zx 52 days, and not well constrained thereafter. An ex- 
trapolation of the cm-band spectrum to the mm-band 
under-predicts the 97.5 GHz flux density by ~ 20- 
5096. Furthermore, the mean spectral index between 
the ALMA (3mm) and SMA (1.3mm) bands over this 
period is fairly flat, Bam = —0.02 + 0.13. Thus, the 


mm-band emission cannot be easily subsumed into this 
additional radio component without either invoking ad- 
ditional high-frequency structure in the emission, or in- 
voking additional sources of systematic uncertainties. 
The spectral index between the cm-band peak in the 
first VLA epoch at ~ 3.46 days at 2.5 GHz and the 
ALMA 3mm (97.5 GHz) observation is 6 = —0.106 + 
0.002, which is shallower than the cm-band spectral 
We test whether phase decorrelation at 


index alone. 


Flux density (mJy) 


TIT TT TT TT TT 


Frequency (Hz) 


Figure 4. The multi-epoch radio SEDs from Figure 5, now 
compared with optical observations and Swift/XRT SEDs at 
similar times, together with a wind FS model (lines). The FS 
cannot explain the observed radio flux and slowly evolving 
radio peak. 


Z 10 GHz in the VLA observations could be responsible 
for a loss in observed flux density at higher frequen- 
cies by self-calibrating the highest-frequency K-band 
(25 GHz) observations in the first epoch at ~ 3.46 days. 
This process significantly reduces the imaging residu- 
als but only marginally increases the 25 GHz flux den- 
sity by 4%, whereas making this flux consistent with 
the 6 ~ —0.1 power law would instead require an in- 
crease of zz 10%. We recognize that all radio observa- 
tions are subject to a systematic flux density uncertainty 
of zz 1096 from the flux calibration process, which is not 
incorporated into the analysis above. Thus, if the true 
cm-band flux were systematically higher by this amount 
at 25 GHz (but not at 2.5 GHz), then there remains a 
possibility that the cm-band and mm-band could yet be 
ascribable to the same additional component. 


We consider the possibility that the entire multi- 
frequency (radio to GeV) afterglow emission may in 
fact arise from a single emission component, but with 
non-standard evolution of break frequencies and fluxes. 
To test this, we anchor the peak of the SED (in F,) 
as observed in the cm band, and evolve it in time ac- 
cording to the inferred peak frequency and peak flux 
evolution from the broken power-law fits to the radio 
SEDs (Table 7). We assume spectral indices of 3, = 2.5 
and $2 = —0.2 below and above the peak, respectively. 
Upon extrapolating this spectrum to the optical, we find 
that an additional extinction of Ay host 7 0.2 mag is 
needed to match the optical flux, although this is still 


13 


only possible at > 3 days. A spectral break is needed 
above the optical in order to not over-predict the X- 
ray flux. We find that an evolution of this break of 
Vp, © 101? Hzx (t/3.46 d)~!3? together with an assumed 
spectral index of 8; ~ —0.86 (to match the NuSTAR 
spectrum) above the break successfully reproduces the 
X-ray light curve at 1keV and 15 keV. An additional 
break is then needed between the hard X-ray and GeV 
bands in order to not over-predict the GeV flux. Fix- 
ing the spectral index above this break to bg ~ —2 
(from the LAT spectrum), we find that the 1 GeV light 
curve can be reproduced by a fixed spectral break at 
Vp, © 20.7 MeV. We plot this model in Figure 6. 

While this fiducial, single-component model ade- 
quately explains the observations in the cm-band, X- 
rays, and at 1 GeV, it under-predicts all optical/NIR 
observations at < 3 days (just like the physical wind- 
model; Figure 2) as well as the entire 1.3 mm (SMA) 
light curve. Furthermore, it does not have any of the 
temporal breaks that are evident in the 3 mm (97.5 GHz; 
ALMA & NOEMA) light curve. It is possible that 
some of these limitations could be resolved by introduc- 
ing additional, potentially moving breaks into the spec- 
tral shape; however, the introduction of these additional 
degrees of freedom would further reduce the predictive 
power of the model and make it even more challenging 
to interpret. Finally, there are no obvious explanations 
for this particular SED shape or evolution of break fre- 
quencies, and thus it is not straightforward to extract 
meaningful physical information from this model at this 
stage. 

In the scenario that the excess cm-band (and possi- 
bly mm-band) emission arises from a separate popu- 
lation of radiating electrons, potential physical sources 
for such a component might be: (i) a radio supernova; 
(ii) a relativistic Maxwellian population of electrons (i.e, 
the “non-accelerated electrons” or “thermal” electrons 
(Eichler & Waxman 2005; Ressler & Laskar 2017; War- 
ren et al. 2018; Margalit & Quataert 2021); (iii) the re- 
verse shock (RS); and (iv) a two-component jet (i.e., two 
FS-like regions of different geometries), possibly with 
energy injection or time-varying microphysical param- 
eters. The radio component is Z 100 times more lu- 
minous than the brightest known radio SNe (Figure 7) 
and peaks significantly earlier (~few days versus ~ 100 
days post-burst); we therefore do not consider emission 
associated with a possible SN to be a viable explana- 
tion for the radio excess in GRB 221009A. The emis- 
sion spectrum from thermal electrons is expected to be 
broad and to cut off steeply above the peak (Ressler & 
Laskar 2017). A detailed test against thermal electron 
models was performed by Laskar et al. (2019a) for the 
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Figure 5. Left: Radio SEDs from our uGMRT, VLA, ALMA, NOEMA, and SMA observations reveal a self-absorbed spectrum 
with a peak at ~ 2 GHz. The spectral index above the peak is shallow (8 ~ —0.2) with a flattening in the mm-band (= 100 GHz), 
as confirmed by extrapolations (open circles) of broken power law model fits (lines) to the observations at < 50 GHz. Right: 
Evolution of peak flux and peak frequency from broken power law fits to the radio SED at < 50 GHz (Table 7) 


Table 7. Radio SED fits and inferred outflow properties 


Time Fpeak Vpeak Bi Ba s R r Emin Beq 

(d) (mJy) (GHz) (1018 cm) (1048 erg) (mG) 
3.46 8640.2 167+0.02 2.5! —0.18 + 0.01 1.6 0.1 1.40+0.03 9.40.1 163+40.02 12.43 + 0.02 
6.44 4.90.8 1.12+0.30 2.5! —0.20* 0.91 £0.27 1434060 6941.4 2.13 40.24 11.9 + 0.05 
12.44 4.30.8 1.19+40.25 2.5!  —0.26 4 0.02 1.8 X 1.2 0.88 +0.12 3.9 40.7 4.02 0.51 18.2 + 0.6 
17.48 2.940.1 0.77+0.06 2.5! —0.23 + 0.01 1.4 X 0.3 1.04+0.11 36+0.2 44.28 0.10 12.2 t 0.1 
28.33 1.90.1 0.56+0.03 2.5! —0.21 Æ 0.01 1.34 0.2 1.00 + 0.08 2.88+0.1 5.30+0.19 10.99 + 0.09 
52.48 1.20.1 0.574012 2.5!  —0.37 + 0.03 0.6 X 0.1 0.58 £0.18 1540.2 7.63 + 0.45 18.3 + 0.6 
76.42 10+0.1 0.45+0.07 2.5!  —0.18-- 0.01 1.07 0.60 +0.14 1.30.2 9.96 + 0.68 15.9 4 0.4 


Nore— Fixed. The last four columns list the equipartition radius, Lorentz factor, minimum energy (in the emitting region), 
and magnetic field, respectively, computed using the formalism of Barniol Duran et al. (2013) (see Section 4.2). 


low-frequency radio excess in GRB 181201A, and they 
confirmed that the radio SEDs in that case were nar- 
rower than thermal electron models would predict. The 
SEDs observed in this case are similarly sharply peaked, 
and thus also unlikely to match our current framework 
of synchrotron radiation from a Maxwellian population 
of electrons. 

In the case of GRB 181201A, the excess radio emis- 
sion was ultimately ascribed to the RS, albeit with non- 
standard parameters. However, the simple power-law 
fits performed above demonstrate that the temporal evo- 
lution of GRB 221009A’s radio component is too slow 
to be ascribed to a RS. To see this, the evolution of 
the spectral peak (va,r) of a Newtonian !' RS in the 


17 The emission can decay slower in a Newtonian RS compared to 


the relativistic RS case. 


ð log Va,r (15g+24)p+32g+40 
Ologt (14g+7)p+56g+28 ^ 


where T x RY is the evolution of the Lorentz fac- 
tor of the post-shock ejecta with radius (Kobayashi & 
Sari 2000). We expect 0.5 « g « 1.5 for a wind en- 
vironment, but even if we set this aside and consider 
arbitrarily large values of g (which results in a slower 
evolution), the expected temporal evolution of the peak 
frequency asymptotes to — ets —0.7 for the value 
of p ~ 1.4 that is required to match the spectral in- 


dex above the peak. Similarly, the evolution of the 


regime Lg, < Var ls 


w 
A 


1—p)/2 
k flux, Fa, = F va OP? This gi 
pea. ux ^, vat ~ vmar | Dmr 1S gives 
Olog Fy a,r g(15p?+42p—2)+4(6p?+14p—5) 
Ologt m 7(29+1)(p+4) 0.94 for 


18 The spectral index in this expression is not 5/2 due to the def- 
inition of Fy,m,r used here, which corresponds to the non-self- 


absorbed flux density at mr. 
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Figure 6. Left: Radio to X-ray SEDs (same as Figure 4), together with a fiducial, single-component, broken power-law SED 


model (lines) where the three break frequencies evolve with time as Vpeak cx 


17039. Vp, X t9 and Uy, CX t°, while the peak flux 


evolves as F, peak x t. ^^ (motivated by the observed radio evolution; Figure 5). The underlying model without the zz 4.3 mag 
of total extinction is shown as thin lines. Right: The corresponding light curves from the GeV to the mm band. The model 
under-predicts the optical observations at < 3 days, the entire 1.3 mm light curve, and does not match any of the temporal 


breaks seen in the 3 mm light curve. 


large g and p ~ 1.4. Both these result in a faster fading 
SED than observed in the cm band for this burst. For 
completeness, we present sample plots of the RS model 
in Appendix B. 

Finally, the radio flux cannot be easily ascribed to FS- 
like emission from a single power-law distribution of elec- 
trons either (e.g., in a two-component jet model). This 
is similar to the scenario explored by Sato et al. 2022 
(although as mentioned earlier, their specific model is 
ruled out by the radio evolution at = 5 d). For px 1.4 
for such a model we find that, using the relations of 
Dai & Cheng (2001) for the regime 1 « p « 2, we 
would expect a peak frequency evolution in this case 
of v, x t™34 for a spherical evolution, and v, œ 109? for 
a jet (i.e., post jet-break), together with peak flux evo- 
lution of Fpeak c t 1? and Fpeak x £71, respectively, 
none of which match the observations. On the other 
hand, while energy injection into the shock producing 
this emission may arrest the decay of this component, 
energy injection by itself would not naturally also pro- 
duce a shallow spectrum above the peak. A combination 
of both p « 2 and energy injection might be a feasible 
match; however, such a model appears somewhat con- 
trived. Whereas time-varying microphysical parameters 
could result in a non-standard evolution of the spec- 
tral peak (usually discussed in terms of vm, although 
the peak is at v4 > Vm here), this would not explain 
the shallow spectrum above the peak either (Panaitescu 
& Kumar 2004; van der Horst et al. 2014; Misra et al. 
2021). 


In summary, to our knowledge, the cm and mm emis- 
sion in GRB 2210094 do not correspond to the evolution 
of any standard emission component, including standard 
prescriptions for forward and reverse shock emission. 
'This is the first time that such a component has been 
captured in such exquisite detail; however, previous ra- 
dio observations of long GRBs have often been sparse, 
leaving the possibility open that similar emission may be 
common in long GRBs. We discuss the potential preva- 
lence of such an emission component in radio afterglows 
of GRBs in Section 4.1. 


4. DISCUSSION 
4.1. Comparison to other GRBs 


We now consider GRB 221009A in the context of 
the broader GRB population. In the X-ray and y-ray 
bands, GRB 2210094 is one of the brightest GRBs ever 
observed (Veres et al. 2022; Ursi et al. 2022; Kennea 
et al. 2022b; Frederiks et al. 2022; Iwakiri et al. 2022; 
Kobayashi et al. 2022; Negro et al. 2023). While this 
is due in part to its proximity, GRB 2210094 is also 
intrinsically among the most luminous known bursts 
at these wavelengths (Figure 7). Our FS modeling 
of GRB 221009A suggests that this superlative lumi- 
nosity is not due to an unusually powerful jet, but 
rather to fortuitous geometric alignment: the beaming- 
corrected jet kinetic energy is average for long GRBs, 
Ex = 50.60+9:99 erg, but the jet’s small opening an- 
gle, Ojet = 1.647928, places it among the most narrowly 
collimated jets in the GRB population (Laskar et al. 
2018). GRB 221009A’s kinetic energy is similar to that 
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Figure 7. Left: Radio spectral luminosity at 8 GHz (rest frame) as a function of rest-frame time for GRBs (deep blue), sub- 
energetic GRBs (light blue), SNe with evidence for relativistic ejecta (dark gray) and hydrogen-poor supernovae (light grey). 
The radio afterglow of GRB 221009A is more luminous than any known radio supernova and comparable in luminosity to the 
afterglows of typical long GRBs. Right: The same plot showing X-ray luminosity in the 0.3-10 keV band. In contrast to the 
radio, GRB 221009Ais one of the most X-ray luminous GRBs for much of its evolution. Adapted from Margutti et al. (2014). 


of other low-redshift GRBs; in particular, it is interme- 
diate between GRB 161219B and GRB 030329 (Figure 
9). Given the high y-ray energy, this implies a relatively 
high prompt efficiency of ny ~ 86%. 

GRB 221009A is also superlative in the quality and 
coverage of radio and mm data obtained. Despite 
GRB 221009A’s extreme brightness at high energies, its 
cm and mm emission is merely average for the GRB 
population (Figures 7 and 8). It is this combination of 
extreme X-ray (and optical) luminosity and mm medi- 
ocrity that makes fitting a FS model to the full dataset 
so challenging: the model overpredicts the mm emis- 
sion at late times unless the mm band remains below 
the peak of the SED, forcing Vm to remain between the 
mm and optical bands throughout the duration of our 
observations and requiring a jet break at ~ 130 days. 

Radio observations can also provide unique insight 
into the physical composition of the jet itself via the de- 
tection of RS emission (Laskar et al. 2013; Perley et al. 
2014; Laskar et al. 2016; Alexander et al. 2017; Laskar 
et al. 2019a,b). However, due to the frequent paucity of 
radio data, deviations of individual data points from this 
basic picture are often ignored, or attributed to other 
effects such as interstellar scintillation (e.g., Alexander 
et al. 2017, 2019; Bright et al. 2019). As discussed in 
Section 3.2 and the Appendix, the exquisite temporal 


and frequency sampling of the radio dataset collected 
for GRB 221009A rule out a standard RS+FS picture, 
as well as other commonly considered model extensions 
(a two-component jet/jet--cocoon, emission from ther- 
mal electrons), with high confidence. We attempt to de- 
termine if similar behavior could have been previously 
overlooked in other well-studied radio GRB afterglows 
in the literature. 

Recently Kangas & Fruchter (2021) systematically 
considered multi-wavelength observations of 21 well- 
studied long GRBs. They showed that while appar- 
ent deviations from a simple FS or FS+RS model in 
the radio within a single event may not appear statis- 
tically significant, when the population is considered as 
a whole, about half of the sample is difficult to explain 
with a standard afterglow model. Levine et al. (2023) 
obtained similar results based on a radio-only analysis 
of a slightly larger sample. In particular, several GRBs 
in Kangas & Fruchter (2021)’s sample (GRB 141121A, 
GRB 160625B, and GRB 171010A) exhibit shallow radio 
SEDs at late times, similar to GRB 221009A. Notably, 
GRB 141121A and GRB 160625B have early radio emis- 
sion consistent with a standard RS model, but the data 
are brighter than the models at late times, similar to the 
issue we faced in attempting to model GRB 221009A’s 
radio emission with a RS. These three GRBs are among 
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Figure 8. Spectral luminosity of GRB 221009A as a func- 
tion of rest-frame time in the 3mm band (red stars), com- 
pared to a population of long-duration GRBs (deep blue), 
supernovae (grey), tidal disruption events (light green), and 
fast blue optical transients (light blue). No K-corrections for 
variations in intrinsic spectral shape have been applied. The 
mm-band emission from GRB 2210094 is a factor of ~ 5-10 
less luminous than GRB 030329 (orange), but of comparable 
luminosity to that of other GRB afterglows. Adapted from 
Eftekhari et al. (2022). 
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Figure 9. Beaming-corrected kinetic energies of GRBs with 
multi-wavelength datasets and known jet opening angles, as a 
function of redshift. GRB 2210094 is both one of the closest 
well-studied GRBs and one of the most narrowly beamed, 
explaining its extraordinary brightness at high energies. 
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the few published events with multi-frequency radio 
coverage extending to such late times, suggesting that 
broad, slowly-evolving radio components like that seen 
in GRB 221009A may be more common than previously 
realized, and may occur in GRBs both with and without 
distinct RS emission at early times. 

While the number of GRBs with well-sampled mm 
light curves remains small, several of them also exhibit 
puzzling behavior relative to other wavebands. For ex- 
ample, GRB 161219B and GRB 181201A’s mm light 
curves were both modeled as the sum of RS and FS emis- 
sion, but fits to the mm light curves in isolation preferred 
a single power law decline with a temporal decay rate 
intermediate between the optical and cm bands, rather 
than a broken power law (Laskar et al. 2018, 2019b). 
Similarly, GRB 221009A’s 3mm light curve evolves as 
~ tbl at late times — steeper than the low-frequency 
radio light curves (~ t^?" at 1.4 GHz), but shallower 
than the late-time optical and X-ray light curves. This 
together with the cm-to-mm SED shape further suggests 
that the mm emission in GRB 221009A may be domi- 
nated by yet another separate component, in addition 
to the FS that dominates the optical and X-rays and 
the slowly evolving component that dominates the cm 
bands. However, standard synchrotron theory struggles 
to produce emission that is narrowly peaked enough to 
dominate the mm while remaining sub-dominant at all 
other frequencies throughout the duration of our obser- 
vations. 


4.2. Equipartition analysis of the radio component 


Assuming that the emission dominating the cm bands 
arises from synchrotron radiation from a shock powered 
by an outflow, we use energy equipartition arguments 
to derive estimates of the emission radius, minimum en- 
ergy, equipartition magnetic field, and bulk Lorentz fac- 
tor of the emitting region. We employ the formalism of 
Barniol Duran et al. (2013), assuming area and volume- 
filling factors of unity, and present our results in Table 7. 
'The emission radius and equipartition magnetic field are 
roughly constant at R z 101? cm and B ~ 13 mG, 
respectively, whereas the bulk Lorentz factor decreases 
from [ z 9.4 to T z 1.3 between 3.46—76.42 days. This 
is lower than the corresponding bulk Lorentz factor of 
the FS (which decreases from T z 73 to T ~ 34 over 
the same period), possibly suggesting some structure 
in the ejecta. The minimum energy increases? from 


19 In the equipartion framework, this apparent increase in Emin is 
due to deceleration resulting in a greater volume of the emitting 
plasma becoming visible, with values computed at later epochs 
providing a stronger constraint on the true minimum energy. 
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Ein & 1.6 x 1048 erg to Emin © 1 x 107? erg over this 
period. 

'The apparent decelerating nature of this component 
yields a lower limit on the Lorentz factor of T'o ~ 9.4, 
which corresponds to an upper limit?" on the ejecta 
mass of Mej 2 ain (T'g6*) y 10-7 Mo for Emin y 
1.6x 108 erg as derived from the radio SED at 3.46 days, 
or Mej S 6 x 1077 Mo as derived from the radio SED 
at 76.42 days. This is similar in magnitude to the in- 
ferred ejecta mass for a typical long-duration GRB with 
beaming-corrected kinetic energy Ex z 5 x 1099 erg 
(Laskar et al. 2014), assuming a typical initial Lorentz 
factor of T'o œ% 10?. This indicates that the outflow 
producing the radio emission in GRB 221009A shares 
characteristics similar to those of standard long-duration 
GRBs. However, the stagnant emission radius inferred 
for this component stands in stark contrast to the 
expanding outflow expected from standard theoretical 
models of a relativistic GRB jet propagating in a smooth 
environment. Finally, the inferred value of the equipar- 
tition magnetic field is higher than that inferred for the 
post-shock B-field in the FS (B ~ 1 mG(t/1 day) ?/4), 
but lower than that inferred for supernova shocks in 
Type Ib/c supernovae (Chevalier & Fransson 2006). 


4.3. Caveats from modeling the FS emission 


We now discuss the mismatch between the FS wind 
model and the data at multiple wavelengths. The FS 
model under-predicts the optical emission at < 1 day 
(Figure 2). Similar excess emission in the past has been 
attributed to RS emission, although the RS in this case 
would not match the radio observations (Appendix B). 
A change in density structure from steep to shallow re- 
mains plausible, although the density profile is already 
steep (k — 2) and the density itself already very low 
(A, z 2 x 1073). Furthermore, such a transition would 
also affect the X-rays, which are in the same spectral 
regime in this model; however, no such transition is ap- 
parent in the light curves. Any additional component 
invoked to explain these optical observations would ei- 
ther need to match the radio SED, or at least, not over- 
predict these. Given the lack of such models, the ob- 
served optical excess is puzzling. 

'The FS model also under-predicts the optical emis- 
sion at Z 30 days. Excess emission above this model 


20 The relationship E ~ TMec? holds only if the system is still 
ballistic. If it has been decelerated by the environment (as the 
stationary radius would indicate), then the Lorentz factor should 
no longer be connected to the ejecta mass, but only to the energy 
and density. Since the Lorentz factor decreases with time, and 
because we infer an increasing energy and decreasing I, the first 
epoch yields the tightest constraint and an upper limit on Mej. 


is also apparent in the Swift/XRT light curve at 1 keV 
in the form of a bump at 30-60 days, and in the final 
NuSTAR epoch relative to the single power-law fit to 
the 15 keV light curve (Figure 1). In fact, a similar ex- 
cess over the broken power-law fits is also apparent in 
the radio light curves on a very similar timescale (Fig- 
ure 1), possibly indicative that the radio emission does 
arise from the same emission region as the X-rays and 
optical, at least in this narrow time interval. The lu- 
minosity of this excess at 11 GHz, the optical 7/-band, 
and at 1 keV is zz 2 x 10?? ergs-!, = 5 x 10*! ergs-!, 
and z 107? ergs-!, respectively. The achromatic na- 
ture of this bump is suggestive of a hydrodynamic effect, 
such as the appearance of a counterjet, a second emis- 
sion component (e.g., cocoon or SN emission), a density 
enhancement (possibly including an encounter with the 
wind termination shock), or energy injection. We ex- 
pect [ zz 40 at ~ 60 days, so the counterjet possibility 
appears unlikely. The luminosity of this component is 
greater than any known SN at each of these wavelengths 
(Figure 7), making such an origin also unlikely. Previous 
work on the impact of density enhancements on optical 
light curves of GRB afterglow suggests that such effects 
lead to minimal deviations in the light curves owing to 
the highly relativistic nature of the jet, even if the en- 
hancement is extreme (Nakar & Granot 2007; Uhm & 
Beloborodov 2007; van Eerten et al. 2009, 2010; Gat 
et al. 2013; Geng et al. 2014). Whereas energy injec- 
tion could conceivably create such a signature, it would 
need to be coupled with a jet-break in order for the light 
curves to not asymptote to a higher flux level. In their 
study of energy injection in GRBs, Laskar et al. (2015) 
find that in 3 out of 4 cases studied, the jet breaks were 
within a factor of ~ 2 from the end of the period of 
energy injection. If the observed bumps in the light 
curve for GRB 2210094 are related to energy injection, 
this would continue such a trend. However, the inability 
of the FS model to directly connect with the radio ob- 
servations makes further investigation of this possibility 
challenging. 

The best-fit wind model over-predicts the observed 
NIR (JH K-band) flux by z 15-30%, with the offset 
worsening with time. On the other hand, no such offset 
is apparent in the optical fits. We find that this is due 
to an apparent change in the spectral slope between the 
NIR and optical bands from SytR-opt = —0.76 + 0.04 
at ~ 4.4 days (Section 3.1) to BNIR-opt = —0.52 + 0.04 
at ~ 17.5 days. Such a hardening of the NIR-to-optical 
spectrum is unexpected and cannot be understood in 
the standard synchrotron framework in our preferred 
spectral regime of Vm S Vopt < vx < Ve. Whereas 
such a color evolution from red to blue is similar to 


that expected from the emergence of an underlying SN, 
inclusion of such a component would over-predict the 
optical observations. Suppressing the afterglow contri- 
bution to both the optical and NIR at this time would 
require either a steepening of the electron index, p with 
time (which would then affect the X-ray fit) or a faster 
evolution of vm, which would worsen the tension in the 
mm-band at ~ 100 days, making this interpretation un- 
tenable. It is also possible that this apparent change in 
ÜNrR-opt instead arises from systematic calibration er- 
rors in the data. Further investigation of this anomaly 
requires better light curve coverage of the NIR afterglow 
and is beyond the scope of this work. 

'The model fit appears slightly worse in the LAT 1 GeV 
band. This is partially due to a slightly shallower decay 
index — the model predicts apar = (2 — 3p)/4 ~ —1.40, 
compared to the observed value of apap £ —1.47. It 
is also in part due to the higher flux in the model (by 
& 50%) compared to the data. The flux offset could 
be addressed by a slightly lower value of v., or, alter- 
natively, both of these could be remedied by a slightly 
larger value of p. Since ve depends on all four physical 
parameters, it is non-trivial to discern the origin of ad- 
ditional constraints on its value imposed by the rest of 
the data. The challenge with a higher value of p, on the 
other hand, lies in fitting the optical-to-X-ray spectral 
index; increasing p while keeping the model flux in the 
X-rays fixed would result in the model over-predicting 
the NIR even more. Thus the current value of p is a 
compromise between strongly over-predicting one or the 
other of the NIR or GeV emission. Finally, we note that 
the LAT spectrum shows some evidence for steepening 
in the last two bins at = 31 ks (= 0.4 days). This could 
be due to the lower photon flux reducing the fluence of 
the highest energy photons in this bin. However, if this 
reduction is real, it could also indicate the movement of 
the maximum synchrotron frequency, Vmax into the LAT 
band. We defer the discussion of this potential effect to 
future work. 

Finally, we speculate on a possible unifying underly- 
ing cause for some of these issues. One possible way 
to resolve the tension between our observations and our 
analytic model might lie in constructing a more realistic 
distribution of relativistic electron energies. Standard 
synchrotron theory assumes that the population of radi- 
ating electrons are accelerated into a simple power-law 
distribution of energies with an abrupt low-energy cutoff 
at Ymin, resulting in a spectral break (vm). Our current 
model for GRB 2210094 requires Vm to remain between 
the mm and the optical throughout our observations; we 
never observe a break in any light curve attributable to 
the transition of Vm through that frequency band, al- 
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though we infer its presence in the optical at S 2 days 
by dint of the constraint on the flux of the v!/? segment 
imposed by the mm-band data at ~ 100 days. Changing 
the low-energy end of the relativistic electron distribu- 
tion would change the SED below Fpeak, exactly the 
regime where our current model struggles to reproduce 
our observations. We defer further exploration of this 
and other possibilities to future work. 


4.4. Synchrotron Self- Compton predictions 


Multiple experiments reported the detection of very 
high energy photons from GRB 221009A (Huang et al. 
2022; Dzhappuev et al. 2022), making GRB 221009A 
the newest member of the very small class of GRBs 
with detected VHE emission (GRBs 180720B, 190114C, 
1908294, and 201216C; MAGIC Collab. 2019; Abdalla 
et al. 2019; Blanch et al. 2020; H. E. S. S. Collaboration 
et al. 2021). GRB 2210094 also exhibited the highest- 
energy photon yet associated with any GRB (18 TeV; 
Huang et al. 2022). Armed with a predictive model, 
however imperfect, we consider whether these VHE pho- 
tons could arise from FS emission. Huang et al. (2022) 
report 5000 VHE photons in the span of ~ 2000 s follow- 
ing the Fermi/GBM trigger, corresponding to a flux of 
z 3 x 1071? ergs^! cm? at 1 TeV (computed assuming 
a LHAASO collecting area of 1 km? and a mean photon 
energy of z 1 TeV, ignoring spectral corrections). At a 
time of ~ 500 s after the burst, we compute both the 
synchrotron spectrum from our FS model, and the corre- 
sponding synchrotron self-Compton (SSC) emission ex- 
pected, the latter by integrating the synchrotron spec- 
trum over an electron distribution back-calculated from 
the locations of the break frequencies at this time. We 
find that a synchrotron flux and SSC flux, at 1 TeV, of 
zx 1.2x 10 9 ergs ! cm7? and = 3 x 1071H ergs"! cm, 
respectively. Thus, the synchrotron spectrum will dom- 
inate over SSC at 1 TeV, and hence the intrinsic spec- 
trum is expected to be ~ —p/2 ~ —1.1 in F, in the 
absence of a high-energy cutoff in the electron spec- 
trum. The synchrotron flux is a factor of ~ 40 higher 
than a rough observed flux computed above; absorp- 
tion due to y-y pair production against the Extragalac- 
tic Background Light (EBL) is expected to attenuate 
the observed spectrum (although, see also Sahu et al. 
2023), and this deficit is of the same order of magnitude 
as, but smaller than, that inferred for GRB 190114C 
(MAGIC Collab. 2019). Thus, it is conceivable that 
the VHE emission for this GRB was produced by the 
FS. We conclude with two caveats. First, this is only 
an order-of-magnitude estimate, and have not consid- 
ered, for instance, a high-energy cutoff in the electron 
spectrum in this calculation. Second, our FS model has 
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several shortcomings, and the true FS model flux at the 
time of the VHE detection remains somewhat uncertain. 


5. CONCLUSIONS 


We have presented multi-wavelength observations of 
the superlative GRB 2210094, spanning fifteen orders of 
magnitude in frequency and four in time. We find that 
GRB 221009A’s NIR, optical, X-ray, and y-ray emission 
can be well-modeled as a synchrotron FS from a highly 
collimated relativistic jet interacting with a low-density 
wind-like medium. GRB 221009A’s high brightness as 
observed from Earth can be attributed to a combination 
of its relative proximity, its large intrinsic luminosity, 
and its high degree of jet collimation combined with an 
on-axis orientation. 

GRB 2210094 strongly demonstrates the need for ad- 
ditional theoretical work to fully understand the ultra- 
relativistic jets seen in long GRBs. While a simple 
FS model is broadly consistent with a large fraction of 
our data, the radio and mm emission in particular are 
difficult to explain within the scope of standard syn- 
chrotron theory. We consider two possibilities: (i) that 
the radio emission is due to an additional synchrotron 
emission component (with a possible second additional 
component required to explain the mm emission) or (ii) 
that our basic analytic models of relativistic synchrotron 
emission need to be modified in some fundamental way. 
We find that the temporal evolution and spectral shape 
of the cm emission are inconsistent with standard ana- 
lytic models for FS or RS emission propagating in con- 
stant density or wind-like media. However, the peak fre- 
quency and peak flux density of this component evolve 
simply with time (as Vpeak c t705 and Fpeay e t^? 
respectively), perhaps suggesting that an analytic de- 
scription of this component might be possible if a non- 
standard assumption is made (e.g., evolving microphysi- 
cal parameters). While fully exploring extensions to the 
standard synchrotron afterglow models is beyond the 
scope of this paper, we briefly speculate that a more re- 
alistic treatment of the low-energy end of the relativistic 
electron distribution may solve some of the issues. 

GRB 221009A's proximity means that it will remain 
detectable with a wide variety of radio facilities for 
years to come, providing a testbed for future theo- 
retical work and an opportunity to further refine the 
synchrotron model applied in this paper. In addition, 
GRB 221009A’s radio brightness and longevity will pro- 
vide rich opportunities for additional science, includ- 
ing directly measuring the physical size of the afterglow 
with VLBI observations (previously only convincingly 
demonstrated for GRB 030329) and constraining the 


magnetic field structure of the jet with polarization ob- 
servations. 
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APPENDIX 


A. AN ISM MODEL 


In this section, we consider an ISM model with p « 2 
in the regime v; < Vopt < vx with an early jet break 
(tiet 0.55 days) in order to match the optical and X- 
ray spectral index and light curves, as discussed in Sec- 
tion 3.1. The low value of p is required to match the ob- 
served NIR-to-X-ray spectral index of Gyrr_-x ~% —0.70 
in the stipulated regime of Ve < Vopt < vx. We set 
p — 1.5 and tune the parameters to match the output 
X-ray and optical light curves. To calculate the light 
curves, we replace € in Granot & Sari (2002) by €e, 
and note that there will need to be a change in elec- 
tron spectrum at some high Lorentz factor in order to 
keep the total energy in accelerated particles finite. In 
this model, ve ~ 3 x 10? Hz < voy, as required in 
order to satisfy BNIR-opt ^ ÜNrR-x (Section 3.1). We 
find that in this model the jet becomes non-relativistic 
at tnr 7 11 days and the resultant model light curves 
over-predict the radio SEDs (Figure 10). This is funda- 
mentally because it is not possible to match the radio 
SED onto the optical with a single synchrotron emission 
component without invoking additional spectral breaks. 
Thus, the ISM, single-jet model is ruled out for this 
burst. 


B. NEWTONIAN RS MODEL 


The excess radio emission described in Section 3.2 
cannot be easily ascribed to RS emission. To demon- 
strate this, we combine the FS model described in Sec- 
tion 3.1 with an RS model with the following param- 
eters: Var & 2.3 X 10? Hz, Vmr & 1.2 x 10H Hz 
and P, ~ 94 mJy, selected to achieve an approx- 
imate match to the cm-band SED at z 17.48 days. 
The RS cooling break is only weakly constrained to 
v... (1day) = 101? Hz in this model, so as to not strongly 
affect the cm-band SED. We also require p ~ 2 in order 
to match the shallow cm-band spectral index above the 
peak at > 2 GHz. We present a Newtonian RS model 
with g — 5 in Figure 11. A higher value of g leads to 
a slower RS evolution. For a wind medium, we expect 
0.5 € g € 1.5. Even with g ~ 5, the evolution of the 
model SED is not slow enough to match the observa- 
tions. A relativistic RS would evolve even faster. Thus, 
the presently available suite of RS models cannot match 
the radio observations for this burst. 
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